Quantum pumping: The charge transported due to a translation of a scatterer 
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The amount of charge which is pushed by a moving scatterer is dQ — —GdX, where dX is the 
displacement of the scatterer. The question is what is G. Does it depend on the transmission go 
of the scatterer? Does the answer depend on whether the system is open (with leads attached to 
reservoirs) or closed? In the latter case: what are the implications of having "quantum chaos" 
and/or coupling to the environment? The answers to these questions illuminate some fundamental 
aspects of the theory of quantum pumping. For the analysis we take a network (graph) as a model 
system, and use the Kubo formula approach. 
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Consider an Aharonov-Bohm ring with a Fermi sea of 
non-interacting spinless electrons as in Fig. lb. Assume 
that the ring is either disordered or chaotic, and that the 
temperature is known. In such a case the ring has a well 
defined Ohmic conductance G"*"". This means that if we 
change the magnetic flux $ through the ring, then we 
have the electro-motive- force — <i> and the current is Z = 
QahrcL ^ (—<!)). Therefore the charge which is transported 
is dQ = — G'°'""d$. Next we can ask what is the charge 
that is transported if we vary some other parameter (X) 
that controls the potential in which the electrons are held 
(in practice it may be a gate voltage). Then, in complete 
analogy with Ohm's law, we expect to get dQ — —GdX, 
where G is a generalized conductance. The calculation of 
G is the so-called problem of "quantum pumping" . 

Most of the studies of quantum pumping were (so far) 
about open systems (Fig.le). Inspired by Landauer who 
pointed out that G"*"" is essentially the transmission of 
the device, Biittiker, Pretre and Thomas (BPT) have de- 
veloped a formula that allows the calculation of G using 
the S matrix of the scattering region 0, 01 • It turns out 
that the non-trivial extension of this approach to closed 
systems involves quite restrictive assumptions 0. Thus 
the case of pumping in closed systems has been left un- 
explored, except to some past works on adiabatic trans- 
port [3, 111 ■ Yet another approach to quantum pumping 
is to use the powerful Kubo formalism [aj l3 ■ In this Let- 
ter we report the first non-trivial demonstration of this 
formalism. Namely, we illuminate the interplay of "quan- 
tum chaos" with non-adiabaticity and environmental ef- 
fects, and in particular we derive specific results for the 
(generalized) conductance in closed (chaotic) system. 

To be specific we ask what is the amount of charge 
which is transported if we make a displacement dX of 
some scatterer or of some wall element ( "piston" ) . The 
answer to this question in case of the ID system of Fig. le 
is well known. Using the BPT formula one gets 
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where kp — mvp/h is the Fermi wave number, and 170 is 
the transmission of the scatterer. This result is analogous 
to the Landauer formula G°'"" = {e'^/h)go. The charge 
transport mechanism which is represented by Eq. has 



a very simple heuristic explanation, which is reflected in 
the term "snow plow dynamics" On the other hand 
in case of Fig. la or Fig.lc, if we translate the scatterer 
at some constant velocity X, it is clear that 
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In Eq. (|2Jl the transmission of the scatterer is assumed to 
be go < 1 (in Fig. la the role of go is played by the relative 
size of the scatterer). Irrespective of go the steady state 
is a distribution that moves with the same velocity as 
the scatterer. (In the moving frame of the scatterer the 
steady state is a standing wave). Though go does not 
influence the value of G, it should be remembered that 
in the limit go ^ 1 it takes an infinite time to get into 
the steady state. 

What about Fig. lb? Here the system has an inher- 
ent time scale Td which governs the relaxation to an er- 
godic distribution in the laboratory frame irrespective of 
the driving. The scatterer should push its way through 
the ergodizing distribution, and therefore its relative size 
(or its transmission) matters. Thus we expect to find 
G oc (1 — go) in the analogous network model of Fig. Id. 
Moreover, we expect to have a dependence of G on the 
overall transmission gx of the ring. Using the Kubo for- 
malism we are going to derive 
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and to discuss the conditions for its validity. In par- 
ticular we are going to discuss whether the effect of 
non-adiabaticity is to induce a crossover from Eq.© to 
Eq.(^. We are going to present both semiclassical and 
quantum mechanical derivations, and to consider the gen- 
eralization of "universal conductance fluctuations" . 

Our model system Q is the network of Fig. Id. It is 
composed of ID wires ("bonds"). The scatterer is rep- 
resented by a delta function which is located on some 
selected bond. Thus the Hamiltonian H incorporates 
a potential U{x) = X{fi^ /2m) S{x — X), where X is 
the location of the scatterer along the bond. The pa- 
rameter A determines the transmission of the scatterer 
go-(l + (A/2fcf)2)-i. 
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FIG. 1: A scatterer (represented by a black circle) is trans- 
lated through a systems that has a Fermi occupation of spine- 
less non-interacting electrons. In (a) the system is a simple 
ring. In (b) it is a chaotic ring (Sinai billiard). In (c) and 
in (d) we have network systems that are of the same type of 
(a) and (b) respectively. In the network, the scatterer ("pis- 
ton") is a delta function (represented as a big circle) located 
a.t X — X. The current is measure at a section (represented 
by a dotted vertical line) through x — xi. In (e) we have an 
open geometry with left and right leads that are attached to 
reservoirs that have the same chemical potential. 

The generalized conductance G is in fact an ele- 
ment of the generalized conductance matrix, that gives 
the response of the current X to the driving via X. 
With X we associate the operator T = —{dH/dX) = 
X{n^/2m)S'{x - X). Since dX is a displacement, it fol- 
lows that is a Newtonian force. We measure the cur- 
rent through another section x = xi along the same bond 
(see dotted line in Fig. Id): 

I ^ — {S{x - xi)p + pS{x - xi)) (4) 
zm 

The Kubo expression for G can be expressed using a gen- 
eralized Fluctuation-dissipation (FD) relation ,7J: 

/•oo 

G = giEp) / CEAr)dT (5) 
Jo 

where Ce{t) is the cross correlation function of X(r) and 
J^(0), and g{E) is the one-particle density of states at 
the Fermi energy. Eq. (jSJ is the Kubo analog of the BPT 
formula. In fact the latter can be obtained as a special 
limit of the former |l7|. Eq.© would become the stan- 
dard FD relation if we were looking for G"*"". In such 
case the symmetric current-current correlation function 
would be involved, and consequently one would obtain 
(^ohm _ lg(^EF)CEp{i~o^Q)- But in the case of our pump- 
ing calculation Ce{t) is antisymmetric, and therefore 
such simplification is not possible. 

Disregarding the assumption of having Fermi occupa- 
tion, Eq. (O has a classical derivation. In such case Ce (t) 
is the classical cross correlation function, and there is an 
implicit thermal averaging over Ep. The analogous quan- 
tum mechanical (QM) derivation is based on perturba- 
tion theory. Later we are going to discuss its limitations. 



In the QM context Ce (t) is the symmetrized correlation 
function. Furthermore, considerations that are based on 
a Markovian treatment of perturbation theory to infinite 
order, suggest the replacement 

Ce{t) ^ CB(T)e-(r/2'')- (6) 

where F cx is related to the non-adiabaticity of the 

driving Q , or more generally it may incorporate also the 
influence of an external bath. Thus we get 

where 

CEH = T7^[C{E+huj,E)+C{E-nuj,E)] (8) 

and 

C{E',E) = 2nY,I,,r,AE'~E,n)T^n5{E-E^) 
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= - trace [I Im[G(£:')] Im[G(i^)]] (9) 
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In the first equality C [E' , E) is expressed using the 
eigen-energies En and the matrix elements of T and T. 
In the second equality it is expressed using Im[G] — 
-ii(G+-G-) where G± = l/{E-n±iQ) are the Green 
functions of the system. 

Thus we have three options for calculation. In the 
classical treatment it is simplest to calculate directly the 
correlation function in the time domain. In the QM case 
we can use expressions for the Green function in order 
to make the calculation. Optionally we can express the 
conductance using the eigen-energies and the matrix el- 
ements of the associated operators: 

where f{E) is the Fermi occupation function. The latter 
expression is valid also in the strict adiabatic limit F — > 
where it can be regarded as geometric magnetism Q. 

We can get a semiclassical estimate for G by study- 
ing the classical correlation function C(t). But first we 
should define what classical calculation means in the con- 
text of this network model. Recall that dX is displace- 
ment, so J- has the meaning of Newtonian force. There- 
fore in the classical calculation J-{t) consists of spikes 
whose area has the meaning of impact (= 2mui?). Pos- 
sibly it is more intuitive to think of the scatterer as a 
rectangular barrier with two vertical walls. The verti- 
cal walls of the scatterer are regularized by giving them 
finite slops. In such case the spikes of J-{t) become 
short rectangular pulses of some duration tq and height 
(2m?;i;')/To. Obviously this regularization drops out from 
the final result, because the product 2{t)F{0) is weighted 
by the probability of having non zero ^(0), which is 
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ivFTo)/{2L). The possibility to tunnel through the scat- 
terer is taken into account by adopting a stochastic point 
of view. Namely, upon collision there is a probability go 
to go through the scatterer (in such case there is no im- 
pact). Using the above stochastic picture one deduces 
that the short time correlations are 



C{r) = ^e{2mv,) 



(11) 



where ri — {xi—X)/vf- However there are tails due to 
multiple reflections, leading (after geometric summation) 
to 



C{T)dT = -e^ 
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Substitution into Eq.lO leads to Eq.Q. The validity of 
the final result can be double checked by solving classical 
Master equation to find the (quasi) steady state solutions 
of the problem. The current in the steady state, to linear 
order in the rate of the driving, leads to the same result 
(Eq.(|3Jl) for the conductance. 

We turn now to the proper quantum mechanical calcu- 
lation. The matrix elements between eigenstates of the 
network are: 

Turn = -«|^ (^-"9^" - aV;"V"),^,^ (13) 
Tnm = -A— (V"9V" + aV^"^'"),-^ (14) 

2m 

where the gradient dip should be interpreted as the av- 
erage value of the left and right slopes. [To derive this 
result it is convenient to regard the delta function as a 
narrow rectangular barrier.] Without loss of generality 
we set from now on X = 0. It is convenient to express 
^nm using the wavefunction at x — +0. Thus we get 



= -A— irdip"^ + - >^r^n.=+o (15) 

Substitution of H13|) and H15|) into (|10|l leads to an expres- 
sion that can be written in terms of the Green function G. 
[An alternate procedure is to substitute in Eq. © the im- 
plied differential representation of the operators] : 



G = 



( - — ] {G^Gi^, — G^^iG^ + G:^G^/ 

TT V 2m 



-G^,Gi + XiG'J,G' -G"Gi,))E (16) 



where = — J dE f'{E) implies thermal averaging 

and we have defined = Y{.eG{x,x']E + iT/2), and 
= Im G(a;, x'; E + iO). The subscripts indicate deriva- 
tives with respect to x and x' . The expression is evalu- 
ated for X = +0 and x' = xi. 

In case of a network the Green function is given by 9] 



G(x,x';£;) = - — 



A, 



(17) 



where the sum extends over all the paths that start at x 
and end at x' . The are the product of the associated 
transmission and reflection amplitudes (t; or for each 
encountered vertex i), while L.p is the total length of the 
path. Upon substitution in l|16() we get a double sum 
^^^^ paths with endpoints x — -fO and x' = Xi. In 
a term that involves derivatives the amplitude Ap (Aq) 
is multiplied by a sign factor s and/or s', which indi- 
cates respectively the initial and final sign of the velocity. 
Gathering all the contributions, one ends up with 
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where Lr — 2hvF/^- Above we neglected off diago- 
nal terms that involve pairs of trajectories with different 
lengths L. This is justified if the energy averaging is over 
a sufficiently large range. It is also important to realize 
that any trajectory that starts at cc = -1-0 and departs in 
the positive direction represents in fact the contribution 
of two degenerate paths: one starts with a positive ve- 
locity, while the other starts with a negative velocity but 
is immediately reflected. Assuming that this is the only 
significant length degeneracy we get 
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The sum is over the same paths as in Ea. (fT7|) . and it can 
be verified that the above mentioned degenerate paths 
adds correctly. The summation over p involves a geomet- 
ric sum in [5t — (1 — 5t)] and gives the factor 177-/(1 — gr)- 
Thus we see that a careful treatment within the frame- 
work of the diagonal approximation recovers the classical 
result. 

We turn to discuss the validity of our result. The 
derivation of the Kubo formula Ea. (|10|) assumes F ^ A(,. 
By definition the bandwidth Ab « fi/Td is the energy 
range \En — Em\ < A;, for which Inm^mn are non neg- 
ligible. It is determined by the classical correlation time 
Tci that characterizes C(t). For a generic chaotic sys- 
tem the mean level spacing is A oc /i"^, where d is the 
dimensionality of the system. Hence the bandwidth in 
dimensionless units is 6 = Af,/A cx; Ti^~'^ . It follows that 
& ;» 1 is the generic case for any quantized chaotic sys- 
tem. For a chaotic network d = \, and b is roughly equal 
to the number of bonds. 

There is a practical implication of the the above discus- 
sion to pumping in general. For some geometries we have 
Tci <S Ti . The notable example is the dot-wire geometry 
of Ref.Q where is related to the motion inside the dot, 
while the current is measured outside at a section on the 
(very long) wire. Thus we may have F > Ti/ti in Eq.@ 
without breaking the validity condition F < Ti/tci- Con- 
sequently we stay only with the short time correlations 
in Ea. lfTT|l . and we get Eq.(P) rather that Eq.©. 

An additional assumption enters into the derivation 
of the generalized FD relation and the associated Green 
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FIG. 2: Comparison of the exact quantum result Ea. ljlOII for 
(G) with the diagonal approximation Eg. 1181 in case of the 
network of Fig. Id. The average is taken over more than 20000 
levels around Ep, while the calculation (for each Fermi level) 
was performed in an interval of 32000 levels. The valency v 
of each vertex is picked up randomly. The transmission of the 
piston is go ~ 0.1. The perpendicular dotted line indicates the 
border of the regime where the QM calculation is valid (see 
text). We also plot the standard deviation SG/{G), while the 
inset displays the distribution P{G-{G)) for F = O.OOOIA. 
Notice the slight asymmetry. The smooth line is the best fit 
with 0.068/(0.024 + \x + 0.065|^'^). In these calculations the 
temperature is T = 0. 

function expression. Namely, it is assumed that F ^ A. 
In order to test the significance of this assumption, we 
consider a generic network with 6 = 35. In Fig. 2 we plot 
the exact result for G as a function of F. We also plot 
the dispersion of G. One may regard this dispersion as 
the analog of universal conductance fluctuations. As F 
becomes larger than A these fluctuations are smoothed 
away. We also see that the result is quite insensitive to 
the exact value of F as long as F ^ A;,, which is the 
regime where the quantum mechanical derivation of the 
Kubo formula makes sense. Throughout all regimes the 
diagonal approximation is very precise, as expected for 
quantities related to the statistics of matrix elements. 

If the condition F <^ At breaks down we enter into a 
non-perturbative regime where the QM recipe Eq. @ does 
not hold. However, if the system has a classical limit, 
then it can be argued (on semiclassical grounds) that in 
the non-perturbative regime the classical calculation can 
be trusted. Since the classical calculation gives the same 
estimate as the diagonal approximation, it follows that 



there should be no apparent breakdown of validity as 
we cross from the perturbative to the non-perturbative 
regime. We have studied this issue in the context of 
energy absorption (the G°*"" conductance). 

It is now interesting to discuss what happens in the 
non generic case 6 ~ 1. For this purpose we consider two 
specific examples: a scattcrer on a closed ring (Fig.lc); 
and a scatterer on a disconnected bond (Fig.le but with- 
out the reservoirs). In the first example the adiabatic 
result Eq.® can be recovered from Eq. © by substi- 
tuting gx — go, while in the second case G = because 
gx = 0. The latter result requires further discussion. 
Since gx = the possibility of getting a steady state cir- 
culating current is blocked. The zero order adiabatic pic- 
ture is as follows: At any moment the lowest energy levels 
are populated, hence at any moment the charge distribu- 
tion is roughly uniform (crgodic). If we plot the current 
as a function of time, we find that the "snow plow" dy- 
namics is counter balanced by adiabatic passages of the 
particles through the moving barrier. The latter manifest 
themselves in the current as short spikes that compensate 
the otherwise steady (snow plow) fiow of current. The 
statistical properties of the current should be regarded as 
the simplest example for the generalized universal con- 
ductance fluctuations that we have discussed previously. 
If the driving is non adiabatic, then the particles do not 
have the time to make adiabatic passages through the 
scatterer, and then the "snow plow" dynamics becomes 
more effective. Thus the non-adiabatic translation of the 
scatterer induces a steady non-zero current in the bond, 
which is associated with accumulation of charge on the 
heading side of the scatterer and depletion in the trailing 
side. Hence the current within the bond becomes (during 
some transient period) of the same order of magnitude 
as in the case of an open system (Fig. le), where the 
reservoirs are assumed to be of an infinite size. 

In summary, the Kubo approach to quantum pump- 
ing allows to explore the crossover from the strictly adia- 
batic "geometric magnetism" regime to the non-adiabatic 
regime. In particular we were able to derive specific re- 
sults for the generalized conductance, using either classi- 
cal stochastic modeling or diagonal approximation, which 
are supported by numerical analysis. 
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